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ABSTRACT 

Using high resolution TV-body simulations wc address the problem of emptiness of 
giant ~ 20/i _1 Mpc-diameter voids found in the distribution of bright galaxies. Are the 
voids filled by dwarf galaxies? Do cosmological models predict too many small dark 
matter haloes inside the voids? Can the problems of cosmological models on small 
scales be addressed by studying the abundance of dwarf galaxies inside voids? We find 
that voids in the distribution of 10 12 Ii~ 1 Mq haloes (expected galactic magnitudes 
~ M*) are almost the same as the voids in 10 11 /i~ 1 M Q haloes. Yet, much smaller 
haloes with masses 1O 9 /i _1 M0 and circular velocities u c i rc ~ 20 km/s readily fill the 
voids: there should be almost 1000 of these haloes in a 20/i _1 Mpc-diameter void. A 
typical void of diameter 20ft. _1 Mpc contains about 50 haloes with v C i IC > 50 km/s. 
The haloes are arranged in a pattern, which looks like a miniature Universe: it has 
the same structural elements as the large-scale structure of the galactic distribution of 
the Universe. There are filaments and voids; larger haloes are at the intersections of 
filaments. The only difference is that all masses are four orders of magnitude smaller. 
There is severe (anti)bias in the distribution of haloes, which depends on halo mass 
and on the distance from the centre of the void. Large haloes are more antibiased 
and have a tendency to form close to void boundaries. The mass function of haloes in 
voids is different from the "normal" mass function. It is much steeper for high masses 
resulting in very few M33-type galaxies (v c i rc ss 100 km/s). We present an analytical 
approximation for the mass function of haloes in voids. 

Key words: Cosmology: theory; dark matter; large-scale structure of Universe 



1 INTRODUCTION 

Already more than two decades ago it became clear that 
large reg ions of the Universe are no t occupied by bright 
galaxies iGregp rv fc Thompson! 119781 : iJoeveer et a 
iKirshner et al.lll98llh Large regions of size ~ (10 — 20)Mpc 
devoid of galaxies can be clearly seen in all present deep 
redshift surveys. The observational discovery was soon fol- 
lowed by the theoretical understanding that voids consti- 
tute a natural ou tcome of stru c ture formation via gravita - 
tional instability JPeebleslll98"l iHoffman fc Shaharnlll982l) . 
Together with clusters, filaments, and superclusters, giant 
voids constitute the large-scale structure of the Universe. In 
spite of the fact that the voids are important for understand- 
ing the observed structure of the galactic distribution, they 
attract much less attention as compared with other elements 
of the large scale structure such as clusters of galaxies. 
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lEinasto et alJ (|W89) were the first to estimate sizes of 
voids in different samples of galaxies with measured red- 
shifts. Voids in the CfA redshift catalogs were studied by 
IVogelev et all Jl994l) . iGhigna et alJ " 



1996) made estimates 
of the void probability function (VPF) for the Perseus-Pisces 
region. V PF was estimated fo r the Las Campanas r edshift 
survey by iMiiller et all fcOOd) . lEl-Ad fc Piranl J1997L l200(fl 
studied voids in the Optical Redshift Survey (ORS) and in 
the IRAS 1.2-Jy survey. They found that large voids with 
radius ~ 20/i~ 1 Mpc occupy about 50% of the volume of the 
Univer se. Void distribution in th e PSCz catalog was stud- 
ied by iPlionis fc Basilakos! J2002f) and by iHovle fc Voeelevl 
(2002) with approximately the same conclusions regarding 
the sizes of voids and the fraction of occupied volu me. For 
more detailed review of observational efforts see iPeeblesI 
teOOll) . It should be noted that in spite of significant ef- 
forts, there remain some crucial unresolved issues regard- 
ing the properties of voids. Voids are defined using bright 
high surface brightness galaxies. This is quite understand- 
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able: finding and measuring redshifts for dwarf or low sur- 
face brightness galaxies is difficult. For exa mple, the abso- 
lute m agnitude limit for the sample used bv lVogelev et all 
il994f) was Mb = —19.3 (as suming the Hubble constant 
h = 0.7). The sample used bv lMuller et alJ (l2000t) was lim- 
ited by Mr = —20. In other words, most of the samples 
are probing voids using galaxies comparable with the Milky 
Way. 

There were several attempts to find dwarf galaxies 
in few individual voids llLindner et al]ll 996l: IPooescu et all 
ll997l:lKuhn et alJll997l iGroein fc GeHerlll999l) . The overall 
conclusion is that faint galaxies do not show a strong ten- 
dency to fill up voids defined by bright galaxies. The limits 
on absolute magnitudes of observed galaxies are better than 
for the large samples, but not overwhelmingly so. Fo r exam- 
ple, one of the voids studied bv iKuhn et al.l 1)19971) was at 
a distance of ~ 3000 km/s, but at that distance the obser- 
vational sample was complete only up to Mb = —18.0. For 
other voids the limit was only Mb = —20.0. The strongest 
arguments th at voids are no t populated by dwarf galaxies 
were given bv lPeebies! (1200 ll) who points out that the dwarf 
galaxies in the ORS catalog follow remarkably close the dis- 
tribution of bright galaxies: there are no indications that 
they fill voids in the distribution of bright galaxies. In this 
case the ORS catalog can "see" galaxies with absolute mag- 
nitudes Mb = —15.5 up to 10 Mpc distance and within this 
distance the voids are clearly empty. A potential problem 
with this statement is that we do not know whether the 
ORS catalog is missing or not low luminosity and low sur- 
face brightness galaxies. At these magnitudes the galaxies 
are likely to have low surface brightnesses. 

To summarize, observations indicate that large voids 
found in the distribution of bright (~ M„) galaxies are 
empty of galaxies, which are two magnitudes below M*. The 
situation at lower limits is not clear. 

Void pheno menon was a targ e t of many theo- 
retical studi es feinasto et al, Il99ll: I Sahni et all 
Ghigna et all Il994t Il996l: iFriedmann fc Piran " 



Arhahi-RiHgoli fc Mullet! I200S 
Benson et alJ l2003t lAntomiccio 



iMathis fc Whit. 
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was st udied bv lEinastoet^dTl^9lF and bv lGlrrgna et alJ 
l)l994h . Our main interest is not the statistics or the 
shapes of the voids. We focus on the issue of emptiness of 
large voids. Are voids empty or are they filled with dark 
matter haloes? What is the structure of the dark matter 
distribution in voids? Does it present a problem for the 
standard cosmological model? These are the questions we 
are trying to address in this paper. Emptiness of voids is of 
additional interest in view of problems of the hierarchical 
models on small scales: the large abundanc e of dark matter 
satellites of Milky Way size haloes ijKlvpin et alJ 1999a; 
iMoore et alj|l999f) and problems with explaining rotation 
curves in central par ts of dwarf and low surface brightness 
(LSB) galaxies (e.g . iMoorel 11994 iFlores fc Primacklll994t 
Ide Blok et aD 1200 j Ivan den Bosch fc S waters! 1200 ll) . Both 
problems are on scales of 10 9 - 10 10 /i" 1 M Q , which can be 
probed by the abundance of dwarf galaxies in voids. 

To s ome deg r ee th e interest for dwarfs in voids is in- 
spired bv lPeeblea B2001» . who claims that the CDM models 
have severe problems: they predict too many dwarfs. While 
Peebles did not try to make any quantitative estimates of 
the number of galaxies or dark matter haloes inside voids of 



large size, the reasoning seems to be simple and straightfor- 
ward. At large redshifts the density in a region, which later 
will become a void, is not much different from the average 
density of the Universe at that redshift. Thus, the fluctua- 
tions grow and haloes collapse. At later times the density 
in the region declines and the fluctuations effectively stop 
growing. The number of collapsed haloes is preserved in the 
comoving volume of the void. This leads to a large number 
of expected haloes and galaxies in the void. For example, 
a void with a density of 1/10 of the average density of the 
Universe is expected to have the number density of galaxies 
roughly one tenth the average density of galaxies in the Uni- 
verse, which gives many galaxi es because of the large void 
volume. IMathis fc White! ( 2002) argue that gravity removes 
haloes from voids. That would reduce the number of haloes 
and galaxies. Our results show that this does not happen 
and thus cannot solve the problem. The simple argument 
of stopping the growth of fluctuations and of subsequent 
dilution of the halo density b y void expansio n must work 
at some scales. In that respect fPeeblej i200ll) is right. The 
only question is what is the scale and what happens on larger 
scales. 

Significant progress in u nderstanding the void struc- 
ture was made rec ently by M athis fe White! 1120021) and 
iBenson et alJ £2003), who used a combination of A-body 
simulations with semi-analytical methods to predict abun- 
dance of galaxies in voids in cosmological models. It was 
found that the voids are empty even of dwarf "galaxies". 
Unfortunately, the mas s res olution in simulation s used by 
IMathis fc White! (1200211 and IBenson et all (I2003T) is low if 
one wants to a ddress the iss u e of d warfs. The best simu- 
lations used bv IBenson et alJ (I2003F) had the particle mass 
1.4 x 10 10 /i _1 Mq, which leads to the minimum halo mass 
of few times 10 h~ Mq, which is not much smaller than 
the mass of our Milky Way galaxy. IMathis fc White! <l2002l) 
had better mass resolution of 3.6 x 10 9 /i _1 M o giving the 
minimum halo mass 3.6 x lO lo /i -1 M0. One of the goals of 
our paper is to extend the limit to much smaller masses to 
find what happens with real dwarfs in large voids. Indeed, 
our mass resolution is almost a hundred times better: we are 
able to detect haloes with mass 10 9 /i _1 M Q . We also develop 
analytical estimates of the mass function of haloes, which 
we test against simulations and then apply to much smaller 
masses. 

One signific a nt ad vanta ge of the app r oach used by 
IMathis fc White! J2002I) and IBenson et alJ J2003h is that 
they were able to estimate luminosities of galaxies hosted by 
dark matter haloes. We do not try to estimate luminosities. 
Instead, our high resolution simulations provide the maxi- 
mum circular velocities of haloes, which gives us a good idea 
of what kind of galaxies we are dealing with. We note that 
in any case the estimates of luminosities in semi-analytical 
models are still very uncertain for dwarf haloes: physics 
of these galaxies is still poorly understood. This problem 
is worsened by low m ass r esolution in the sim ulations of 
IMathis fc White! (I2002T) and IBenson et alJ ^2003), who used 
haloes with as few as 10 particles to track the history of 
smallest galaxies. The estimates of the maximum circular 
velocities are better because they do not depend on what is 
assumed about the star formation in dwarf galaxies. Yet ac- 
curate estimates of circular velocities require high resolution 
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TV-body simulations: internal structure of haloes should be 
resolved. 



In order to get a rough estimate of what luminosities 
may be expected for galaxies hosted by haloes in our sim- 
ulations we present examples of dwarf galaxies in the Lo- 
cal Group, which have measured circular velocities and lu- 
minosities. If haloes in simulations host the same type of 
galaxies, we should expect the same luminosities. NGC 6822 
and NGC 3109 are irregulars with circular velocities about 
60km/s and a bsolute magnitude s Mb = —15.8 in the case 
of NGC 6822 jHodee et alJll99lH\Veldrake et aljkoo^ and 
M B = -15.2 for NGC 3109 l)Mateolll99Sft . For haloes with 
the maximum circular velocity of 60 km/s the virial mass 
is about (2.5 — 4) x 10 10 /i _1 Mq if we assume halo concen- 
trati ons in the range 1 — 2 0. For galaxies with this virial 
mass iMathis fc White! (|2002l) give slightly larger luminosity 
of M B = -16.3. Haloes with virial mass (2-3) x 10 11 /i" 1 M a 
play important role for voids. We find that voids start to 
fill up with haloes wi th masses smaller than this mass. 
IMathis fc White! i2002ft give absolute magnitudes Mb be- 
tween — 17.5 and —18.5 for galaxies with these virial masses. 
The maximum circular velocities for haloes of this mass are 
about 100 km/s, which is comparable with the circular veloc- 
ity of the spiral galaxy M 33 (Mb = —18.5), « c irc = 120 km/s 
JCorbelli fc Saluccill200(i|) . 



We investigate the formation of voids in the stan- 
dard cosmological model: a spatially flat A-dominated Uni- 
verse with scale-invariant adiabatic Gaussian fluctuations. 
We use the following cosmological parameters: J7m = 0.3, 
JIa = 0.7, ug = 0.9, and h — 0.7. The age of the Universe 
in this model is approximately 13.5 Gyrs. These parame- 
ters are favored by recent cosmological observations (e.g . 
iFreedman et al l200lURiess et al J200ll:ls^e"r g el et al.l200.-tl. 
The normalization < js = 0-9 of our simulations is a rather 
conservative value jBunn fc White! Il997l : IViana fc Lid"dlel 
1996). Some recent observational results suggest a sub- 
stantially lower normalization o r a lower density param- 
eter Qm jReiprich fc Bohringerl 120021; IViana et alj l2002t 
iBahcall et al.ll2003ft . lPierpaoli et alJ (120031) found <t« = 0.8. 
their Table 1 contains a compilation of recent estimates of 

The plan of the paper is as follows. In Section we 
present our numerical simulations. We discuss briefly the 
void finding algorithm which we used to find voids in the 
distribution of dark matter haloes of a low resolution simu- 
lation and the resimulation of the regions of selected voids 
with a higher mass resolution. Section [3] is devoted to the 
formalism of analytical predictions of the mass function in 
voids. In Section^Jwe discuss the mass and halo distribution 
in voids. In Section we compare the analytical predictions 
of the mass function with the mass function measured in 
the high resolution simulations of voids. We discuss how the 
mass function depends on the normalization and a chang- 
ing slope of the power sp ectrum as recently proposed by 
the WMAP collaboration JSpergel et al] 12003ft . Finally, in 
Section |S| we summarize our results. 



2 NUMERICAL MODELS 
2.1 The code 

The Adaptive Refine ment Tree (ART) A-body code of 
iKravtsov et alj il997f) was used to run all numerical sim- 
ulations analyzed in this paper. This code uses Adaptive 
Mesh Refinement technique to achieve high resolution in 
the regions of interest. The computational box is covered 
with a uniform grid which defines the lowest (zeroth) level 
of resolution. The code then reaches high force resolution 
by recursively refining all high density regions using an au- 
tomated refinement algorithm. The shape of the refinement 
mesh can thus effectively match the geometry of the region 
of interest. This algorithm is well suited for simulations of 
a selected region in a large computational box, as in the 
simulations presented below. During the integration, spatial 
refinement is accompanied by temporal refinement. Namely, 
each level of refinement, I, is integrated with its own time 
step Aai — Aao/2 l , where Aao is the global time step of the 
zeroth refinement level. In addition to spatial and tempo- 
ral refinement, simulations described below also use a non- 
adaptive mass refinement algorithm to increase the mass 
(and c orrespondingly the force) resolution inside a specific 
region jKlvpin et al.ll200lft . 

We start with running a low resolution simulation with 
128 3 particles covering the whole computational box. For 
that we make a realization of the initial spectrum of pertur- 
bations with 1024 3 particles in the simulation box. Initial 
coordinates and velocities of the particles are then calcu- 
lated using all waves ranging from the fundamental mode 
k = 2tt/L to the Nyquist frequency k = 2?r/L x A 1/3 /2, 
where L is the box size and N — 1024 3 is the number of par- 
ticles in the simulation. Then we replace every 8 3 = 512 par- 
ticles with particles of larger mass. A large-mass (merged) 
particle is assigned a velocity and displacement equal to the 
average velocity and displacement of the small-mass parti- 
cles. Once a simulation with 128 3 particles is completed, we 
select voids and identify all particles inside the voids. 

Then we restart the simulation keeping high mass reso- 
lution for particles in voids. Particles outside the high resolu- 
tion region are merged in several steps so that the high reso- 
lution region corresponding to 1024 3 particles is surrounded 
by shells with resolution corresponding to 512 3 and 256 3 par- 
ticles. The remaining part of the simulation box is simulated 
in low mass resolution corresponding to 128 3 particles. High 
force resolution is only achieved inside the region with high 
mass resoluti on. The details of th is multi-mass technique are 
described bv lKlvpfn et al] fcOOll) . 



2.2 Simulations 

In order to study the formation of large voids, the simulation 
box should be sufficiently large; we use a cube of 80/i _1 Mpc 
on a side. This is sufficient because we are not interested 
in a statistics of large voids, which would require a signifi- 
cantly larger volume. Our main interest is in the structure 
of a typical ~ 10/i~ 1 Mpc-radius void. The 80/i _1 Mpc box 
is large enough for that. We focus on the formation of small 
structural elements (haloes and filaments) inside voids, for 
which we need the highest possible mass resolution. 

The limitation of 1024 3 particles gives the 4.0 x 
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1O 7 /i _1 M0 per particle. This leads to the minimum halo 
mass of 10 9 /i -1 Mq. A halo with this mass has the maxi- 
mum circular velocity ~ 20 km/s. Inside the high resolution 
region we reach the force resolution (one cell) of 0.6/i -1 kpc. 
We make 250000 time-steps on the highest resolution level. 

The identification of haloes is always a challenge. The 
widely used halo-finding algorithms, the friends-of-friends 
(FOF) and the spherical overdensity, both discard "haloes 
inside haloes", i.e. satellite haloes located within the virial 
radius of larger haloes. We have developed two algorithms 
that do not suffer from this drawback: the hierarchical 
friends-of-friends (HFO F) and the bound density maxima 
algorithms (BDM, see lKlypin et al.lll999bh . 

The algorithms are complementary. They find essen- 
tially the same haloes. Thus we believe that the algorithms 
are stable and capable of identifying all dark matter haloes 
in our simulations. The advantage of the HFOF algorithm is 
that it can handle haloes of arbitrary shape, not just spher- 
ical haloes. The advantage of the BDM algorithm is that 
it describes the physical properties of the haloes better by 
identifying and removing unbound particles. In particular it 
estimates not only the mass of a halo, but also its maximum 
"circular velocity", u c i rc = \/ GM/R. This is the quantity 
which is more meaningful observationally. Numerically, v c i IC 
can be measured more easily and more accurately than the 
mass. In order to compare the velocity function of haloes 
measured in simulations with analytical predictions one has 
to convert the virial mas ses into circular velocit ies assuming 
an NFW density profile dGottlober et alJll999l) . 

2.3 Finding voids 

In order to identify voids, we start with construction of the 
minimal spanning tree for selected haloes. Typically we se- 
lect haloes with mass larger than 2 x 10 11 /i _1 Mq, but dif- 
ferent criteria were also used. Then we search on a grid 
with mesh size 0.6h~ Mpc for the point in the simulation 
box which has the largest distance Ri to the set of haloes. 
This is the centre of the largest void the radius of which 
is Ri. We exclude this void and search again for a point 
with the largest distance to the set. This gives the second 
lar gest void and so on. The algorithm is s imilar to that used 
bv lEinasto et alJ <1989ME1-Ad fc Piranl Jl997l) use a some- 
what more complicated search algorithm based on "wall" 
and "field" galaxies, where field galaxies are allowed to be 
also in voids. 

In principle, o ur algorithm (as the algorithm of 
lEl-Ad fc Piranl Il997l) allows for the construction of voids 
with arbitrary shape: the starting point is a spherical void 
which can be extended by spheres of lower radius which grow 
from the surface of the void into all possible directions. How- 
ever, in the following analysis we have restricted ourselves 
to spherical voids to avoid ambiguities of the definition of 
allowed deviations from spherical shape. 



3 ANALYTICAL APPROXIMATIONS FOR 
THE MASS FUNCTION IN A LOW 
DENSITY REGION 

In this section we provide a formalism for predicting the 
mass function of dark matter haloes in voids. There are 



different ways of constructing the mass function in voids . 
Usi ng the barrier-crossin g formalism of |f k)nd et all Jl99lf) 
and lLacev fc Coll Jl993t) . iMo fc White! Jl996t) generalized 
the Press-Schechter approximation so that it can be ap- 
plied to over- and under-dense regions. The validity of this 
approach has been ver i fied a gainst TV-body simulations by 
iLemson fc Kauf fmannl Jl999l). A b etter approximation is 
provided bv lSheth fc Tormenl fcOOl ST). We compare pre- 
dictions based on this approximation with our TV-body re- 
sults and find that it does not provide an accurate fit to the 
results of simulations. This motivates us to develop our own 
approximation. 

For completeness, we start with presenting the con- 
strained ST approximation. If p vo id is the mean density of 
matter in a void and a(R) is the rms density fluctuation at 
comoving scale R, then the number density n Cj sT of haloes 
with mass M(R) = AnpbR 3 /3 is given by 



:,ST(M) = 



(I) 



exp 



1/2 gvoid |?VVo 2 )| <rda 
M (a 2 - 0-2)3/2 dM 

[B(a 2 ) - B(a 2 )] 2 ' 



2(o-2 - erg) 



where 



2» 2 ) 



71 = 



(a 2 - a 2 )" d"[B(* 2 ) - B(a 2 )] 
n\ d(a 2 ) n 



(1) 



(2) 



(note the correction of the typo with respect to 
ISheth fc Tormenl 120021 : R. Sheth, private communication). 
Here 



B(a 2 ) = a 1/2 5 c [l + f3(a8 2 c /a 2 )- a ], 
B(a 2 ) = a 1/2 5 [l + P(a5 2 /a 2 )- a }. 



(3) 
(4) 



The parameters a = 0.615, f3 — 0.485 are given by the ellip- 
soidal dynamics, while a — 0.707 is adjusted by comparison 
with simulations. The parameter So is the linear underden- 
sity of the void corresponding to the actual nonlinear un- 
derdensity at the moment at which we measure the mass 
function (z = in o ur case). It is calculated from the spher- 
ical top-hat model JSheth fc Tormenll2002T) : 



1.68647 
1.12431 



1.68647 



1.35 



+ 



(l + <5) 2 / 3 
0.78785 



(•») 



Here S denotes the mean density contrast in the void, 
5 = £Void/£>b - 1 = A(i? void )/riM - 1, where A = 
3M{R)/(4ttR 3 q c ) is the mean density in a sphere of radius 
R. We will assume that the density in a void does not depend 
on the distance from the centre of the void. Later we will 
see that this not exactly true, but it is a reasonable starting 
point. In the expression above, Ob is the background density, 
5 C is the characteristic density for collapse as predicted by 
linear theory according to the spherical collapse mode l (for 
our cosmology, S c = 1.676, see lLokas fc Hoffmanll200ll) . 

The parameter 00 is defined as the linear rms fluctua- 
tion on the scale of the void. It is estimated using the linear 
power spectrum and the top-hat filter with radius Ro de- 
fined by Rq = (l + 5)R 3 , oid , where 7? vo id and 5 are the radius 
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and the mean density contrast of the void. The parameter 
(To is essential for this approach. It provides truncation for 
objects with large mass: as the amplitude of perturbation a 
approaches the amplitude of perturbation of the void, the 
number density of haloes goes to zero. The effect of trunca- 
tion is "felt" even for much smaller objects. 

The formalism of the constrained mass function in the 
version proposed bv lSheth fc TormerJ {2002) is rather com- 
plicated and arbitrary in a sense that the series in equation 
is not well motivated and the parameters of the model are 
adjusted by comparison with iV-body simulations. We intro- 
duce an alternative and in our opinion much more straight- 
forward and natural method of predicting the mass function 
in voids. As we will show in Section |S] it also reproduces 
our simulated mass functions more accurately. Our method 
is based on the assumption that the evolution of matter 
distribution in a void proceeds effectively as it would in a 
Universe with cosmological parameters similar to those of 
the void. We treat the void as a Universe with a density pa- 
rameter fiM.void = A(i? void ) with A(i? voi d) measured from 
the simulations. 

The growth of perturbations in such a Universe is also 
slower than in the whole Universe. In order to take this into 
account we change the normalization of the power spectrum 
in the void by assuming a new value of crg, vo id, related to 
the background erg by 




v c]rc > 1 20 km/s 
v a „ > 1 50 km/s 
v a „ c > 200 km/s 



0.00 



5 10 15 20 

number of voids 



Figure 1. Cumulative fraction of the volume occupied by 20 
largest voids in the distribution of haloes with different limit- 
ing circular velocities. More massive haloes define slightly larger 
voids. 



O"8,void — Cs 



-DQ0 £>void(a = l) 
D(a = l) Dvoid(ai) 



(6) 



where D(a) and D vo id(a) are the linear growth factors of 
density perturbations in the background and in the void 
respectively, normalized so that for Dm = 1 and Qa ~ 
we have D(a) = a. We assume that at some initial en = 
l/(l + jzi) (zi ~ 1000) the rms fluctuations in the background 
and in the void were equal. 

For a flat model with Hm + Da = 1 describin g our back- 
groun d Universe the growth factor is given iSilveira fc Waeal 
1994, corrected for typos) by 



-Dflat (a) 



a 2 



Fi 



1 l 11 j3 "m - 
3 ' ' 6 ' Qm 



(7) 



where 2-F1 is a hypergeometric function. For non-flat cos- 
mologies, such as those of t he voids, D(a) has to be 
calcu lated numerically using (iHeathl Il977t ICarroll et ail 
119921) D{a) = [5n M /(2o/(o))] J a / 3 (a)da, where /(a) = 

[l + fl M (\ - l) + H A (a 2 - l)] ~ 1/2 . We then use uncon- 
strained ST approximation JS'heth fc Tormenlll999ll to find 
the mass function of haloes in voids: 



n ST (M) 



1/2 



A 



1 + 



g b S c da 
X M^dM 6XP 



a 2 

aSc 
~2cr 2 



,1/2 



(8) 



where A = 0.322, p = 0.3 and a = 0.707. Here we take 8 C = 
1.62, which is appropriate for the 'open Universe' parameters 
of the voids i n our simulations Dm. vo id = 0.03 — 0.05 and 
fi A = 0.7 fsee lLokas fc HoffmanlEooih . 



4 VOIDS IN SIMULATIONS 
4.1 Mass distribution in voids 

In the simulation box of size 80 ft _1 Mpc we found 1348 dark 
matter haloes with circular velocities larger than 200 km/s 
and 2518 with circular velocities larger than 120 km/s. We 
search for spherical non-overlapping voids in the distribution 
of the dark matter haloes. The 20 largest voids have radii 
larger than 8.81 /i _1 Mpc when the limiting circular velocity 
of haloes is set to v c i IC > 200 km/s while larger than 8.26 
h~ Mpc for haloes with v c i IC > 120 km/s. In Figure we 
show the cumulative fraction of volume occupied by the 20 
largest voids. The voids in the distribution of more massive 
haloes tend to be larger, but the difference is small as long 
as the difference i n mass is not very big. A similar behaviour 
has been found bv lArbabi-Bidgoli fc Miillerl J2002il . This in- 
dicates that there are only rare cases when a given void is 
divided into two parts if the threshold for the circular ve- 
locity (or mass) of the objects defining the void is reduced. 
This statement is supported by our conclusion that there is 
a tendency to find more massive haloes in the outer part of 
voids. 

Five voids were resimulated with mass resolution 4.0 x 
10 7 h~ x Mq. In the low resolution simulation their radii were 
Rvoid = ll.e/i^Mpc, lO.S/i^Mpc, 9.4/i _1 Mpc, 9.1/i _1 Mpc, 
9.1/i _1 Mpc. With the high resolution mentioned above we 
resimulated regions with 10% larger radii, so that the ob- 
jects which define the borders of the voids have been also 
resimulated. The void finding algorithm assumes point-like 
objects. However, these objects have a certain size. More- 
over, they could be surrounded by satellites with smaller 
masses. Since we do not want to include in our analysis the 
objects themselves or their satellites (which in the sense of 
the definition do not belong to the void) we have assumed 
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Figure 2. The distribution of dark matter inside a large void. 
Bottom panel shows the distribution of matter in a void with 
radius 10/i _1 Mpc at redshift 2 = 0. Top: all particles belonging 
to the void are traced back to redshift 2 = 2 and are shown in 
comoving coordinates. The distribution looks like normal large- 
scale structure in the distribution of galaxies: the haloes are in 
filaments, there are "voids" and most massive objects are at inter- 
sections of large filaments. Yet, the whole region is a void: there is 
no Milky Way size halo in the picture. Numerous knots along the 
filaments are haloes with masses 10 9 /i — 1 Mq. Large haloes have 
masses of few times 10 10 h~ 1 M.q. 

a void radius of 10 ft _1 Mpc or 8 ft _1 Mpc for the large and 
small voids respectively. 

In the bottom part of Figure we show a sphere of ra- 
dius 10 /i _1 Mpc centered on the void of radius 10.8 /i _1 Mpc. 
This void does not contain any halo of mass greater than 
2.0 x 10 11 h' 1 M Q . The progenitor of this void at redshift 
z = 2 is not spherical. It is much smaller in comoving co- 
ordinates, i.e. the density contrast of the void with respect 
to the mean density was much smaller at high redshifts. 
Note, however, that this does not prove the statement that 
voids become more spherical during evolution, because the 
selected voids are spherical by definition at z = 0. The void 
and its progenitor are shown in the same projection in co- 
moving coordinates. The obvious shift of the centre of the 
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Figure 3. The mean density in spheres of radius R centered 
on the centre of each of the five voids. The two thick solid lines 
denote the voids with radius 10/i — x Mpc, the three thin solid lines 
denote voids with radius 8/i -1 Mpc. 

void means that the void not only expands anisotropically 
with respect to the background but also that the whole void 
moves with respect to the comoving coordinate system. It 
is interesting to see that there is already a huge number of 
dense filaments at redshift z — 2 which should be observable 
in the Lyman a forest. 

Inside the void at redshift z = 0we find almost the same 
structures as seen in simulations of large parts of the Uni- 
verse: empty regions, filaments and matter concentrations 
at the points where filaments join. However, all masses are 
scaled down by a factor of several orders of magnitude. In 
the crossing points of filaments we find instead of huge dark 
matter haloes hosting clusters of galaxies only small haloes 
which might host dwarf galaxies. Along the filaments even 
smaller haloes are situated. The three-dimensional distri- 
bution of matter seems to be non-uniform. Large nodes of 
filaments seem to be situated nearer to the border of the 
void than to its centre. This visual impression is supported 
by the distribution of dark matter (cf. Figure |3J and haloes 
(cf. Figure 2J. 

Let us first consider the mean density in spheres cen- 
tered at the void's centre. We express the mean density in 
units of the critical density, i.e. in a way similar to the (con- 
stant) Q parameter: A = 3M (7?)/(47r7i 3 £>c) where M(R) is 
the total mass inside a sphere of radius R centered at the 
centre of the void and q c is the critical density. As one can see 
in Figure [3] inside a void the density increases slightly with 
radius and is typically a factor of 10 smaller than the mean 
density £Im = 0.3. Our voids have sm a ller d ensities than 
those described bv iFriedmann fc Piranl J200lft . who found 
that voids typically have half the mean density. 

Large differences can be seen in the environment of 
voids. Most of the voids are situated in regions of low den- 
sity. Up to the radius of 30 /i _1 Mpc the mean density in the 
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spheres is still well below the mean density of the Universe. 
However, the most prominent structure of the simulation, 
a galaxy cluster of mass 2 x 10 15 h~ l Mq, is close to void 
4. Therefore, the density outside this void rapidly increases 
so that the mean density of a sphere of radius 20 h~ Mpc 
centered on this void is already well above £7m = 0.3. On 
the contrary, the most prominent clusters in the box are at 
distances of 51 /i _1 Mpc resp. 47 /i _1 Mpc from the center of 
void 5. Therefore, the mean density in spheres centered on 
the centre of void 5 remains below 0.2 up to radii 30 /i _1 Mpc 
and reaches 0.3 only for radii above 50 /i _1 Mpc . We did not 
find any significant differences in the inner structure of the 
voids so we conclude that the environment of the voids has 
no influence on the void itself. 

4.2 Haloes in voids 

Visual inspection of the void simulations shows plenty of 
haloes in the void. In Section|2]we have described how haloes 
can be identified in dark matter simulations. At first glance 
we can see that haloes are not homogeneously distributed 
in the void (Figure 0, large haloes seem to be concentrated 
in the outer parts of the void. We want to quantify this 
statement. 

We divided our sample of haloes in the five voids into 
two subsamples containing 164 haloes with circular veloci- 
ties 55 km/s < u c i rc < 120 km/s and 207 haloes with circular 
velocities 20 km/s < u c i rc < 55 km/s. Each void was divided 
into five shells of equal volume. In Figure [I] we show the 
mean number density of haloes in the five shells the radius 
of which is normalized to the void radius. The thick line in 
the Figure corresponds to the haloes with higher circular 
velocities, while the thin one to those with lower velocities. 
One can clearly see that there is a tendency for the haloes 
to concentrate more at the outer parts of the voids. This 
tendency is more pronounced for the more massive haloes 
(thick line) , their density in the outer shell is almost a factor 
of three higher than in the central sphere. 

Down to the limit of 1O 9 /i -1 M0 we find thousands 
of haloes in the simulation. Now we are interested in the 
mass function of haloes in different voids and its depen- 
dence on the mean density in the void. We select all 
haloes inside the assumed void radius -Rvoid = 10/i _1 Mpc 
or -Rvoid = 8/i -1 Mpc depending on the void and estimate 
the mass function of haloes in each of the voids. Figure Q3 
shows the five mass functions measured in the five simulated 
voids. The mean dark matter density in the voids (see Fig- 
ure is about A(7? vo id) = 0.03 for the larger voids with 
Rvoid = 10/i _1 Mpc and A(i? void ) = 0.04 for the smaller 
voids with -Rvoid = 8/i _1 Mpc. 

The different mean (under)densities of the simulated 
voids result in different mass functions. The higher the den- 
sity in voids the higher is also the number density of haloes. 
The number density of haloes in voids is about an order 
of magnitude smaller than in the whole box as expected 
due to mean density in v oids which is also ab out an or- 
der of magnitude smaller dGottlober et alJ l2002). This can 
be also seen in Figure |S| where we have shown the mass 
function of field haloes obtained by the Sheth-Tormen for- 
malism. This mass function is in excellent agreement with 
the mass function me asured in the whole 80 ft - 1 Mpc box 
iGottlober et al.ll2002T) . In Figure|S|we have scaled the mass 
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Figure 4. Mean number density of haloes with circular velocities 
55 km/s < v c - lIC < 120 km/s (thick line) and 20 km/s < n c ; rc < 55 
km/s (thin line) in shells of equal volume (Vv id/5) for five voids. 
The radius of the shell is normalized to the void radius. 

functions in voids with the mean density contrast measured 
in the void with respect to the mean density of matter in the 
universe, A(R vo id) / 'flu ■ Now the scatter between the mass 
functions of the different voids is smaller. Note, that the 
shape of the mass function in voids is steeper than in the 
whole box. One can clearly see that after rescaling the over- 
all mass function with the density in voids, the less massive 
halos are as abundant in voids as the general mass function 
predicts. The more massive halos, on the other hand, are 
deficient in the voids. 



5 THE MASS FUNCTION IN VOIDS 

5.1 Comparison of simulated and analytical mass 
functions 

The top panel of Figure Q shows predictions obtained us- 
ing constrained Sheth-Tormen mass function equations Q- 
@ with A(R void ) = 0.05, 0.04, 0.03 and assuming the size 
of the void R vo id = 10/i -1 Mpc (smaller sizes of voids re- 
sult in steeper mass functions). The analytical results p re- 
sented here differ from those of iGottlober et all J2002T) in 
that they assumed a = 0.5 in equations to fit the 

simulated mass function better while here we keep the value 
a = 0.707 advertised by Sheth and Tormen everywhere. The 
constrained mass functions do not provide good fits to the 
iV-body results. They are too steep and are a factor of 2- 
10 below the simulations. They are especially bad for voids 
with very low densities. 

In the bottom panel in Figure [7| we compare the simu- 
lated mass functions with the predictions based on rescaled 
mass function given by equation JHJ . To make the analytical 
predictions we use the same values of the mean density as in 
the simulated voids: SlM,void = A(R vo id) = 0.05,0.04,0.03. 
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Figure 5. Mass functions of haloes in five simulated voids. The 
two thick solid lines correspond to voids with radius -R V oid = 
10h~ 1 M.pc and mean density A(i? vo ; d ) = 0.03, while the three 
thin solid lines denote voids with radius i? vo id = 8/i Mpc and 
mean density A(iJ vo j<j) = 0.04 (see Figure 13*1. 

For voids with these densities we get the corrected val- 
ues for the normalization of the power spectrum <rg. vo id = 
0.25, 0.20, 0.14 respectively. It is clear that the rescaled mass 
functions provide much better fits to the simulations. Yet, 
the fits are not very accurate. One may try to improve the 
fits by changing parameters A, p and o in equation ||SJ|. After 
all, values of these parameters to some extent were tuned to 
fit iV-body simulations. However, this is not a good way to 
improve the approximation: predictions are q uite stable. For 
examp le, the values of A, p and a suggested bv lJenkins et alJ 
( 2001) do not change our results significantly and do not 
produce better agreement with simulations. 

The low accuracy of the approximation can be traced 
to the fact that we assume a constant density of the dark 
matter inside a void. At the same time simulated voids have 
density, which visibly increases close to the void boundaries. 
The number density of haloes (especially massive ones) is 
very sensitive to the average dark matter density as clearly 
illustrated by Figure |1] This inconsistency in treatment of 
the voids is the main reason for poor quality of the fits. The 
tendency of more massive haloes to concentrate in the outer 
part of the voids is also in agreement with the slowly increas- 
ing cumulative volume fraction of voids defined by samples 
of haloes with different minimum circular velocities (cf. Fig- 
urc-0. In fact, if those haloes were uniformly distributed one 
would more often expect that one of the big voids will be 
divided into two smaller ones if one decreases the threshold 
of the halo mass or circular velocity defining the voids. This 
is not the case as one can see from the small differences in 
the volumes of the largest voids defined in the set of haloes 
with different circular velocities (see Figure 0. 

In order to reduce the effect of the varying density, we 
find the simulated halo mass function using void radii which 



10.00 




Figure 6. The same mass functions as in Figure|2]but scaled by 
A(-R V oid)/^M- The dashed line denotes the mass function of the 
field according to the Sheth-Tormen formalism, thin and thick 
solid lines correspond to voids with different radii as in Fig. [F] 

are 20% smaller than their actual radii. Figure |H| shows the 
results. The full curves represent the mass functions mea- 
sured in the five voids with the reduced radii. They are now 
steeper at the high mass end. The dashed curves show the 
predictions based on the rescaled Sheth-Tormen approxima- 
tion (the same as in the bottom panel of Figure 0. The 
agreement between the analytical predictions and numeri- 
cal results is now significantly improved. 

5.2 Dependence of the halo mass function in 
voids on different parameters 

Even now the normalization of the power spectrum of den- 
sity perturbations has some uncertainties. Observational 
constraints before recent WMAP res ults give uncertainty of 
about 10% for the parameter erg (see lBahcall et ai"ll2003l for 
SDSS estimates and for a discussion of other measurements 
which have a tendency t o predict low values of a$). How- 
ever, the WMAP results jBennett et al.ll2003ft again favour 
as = 0.9 if the perturbation spectrum is scale invariant with 
n = 1. The mass function in voids is affected by all these 
uncertainties: it becomes steeper with decreasing normaliza- 
tions of the background power spectrum (with other cosmo- 
logical parameters unchanged). 

The combination of WMAP results with other CMB 
measurements, 2dF and Lyman a forest results indicate a 
lower normalization of as = 0.84 together with a scale- 
dependent slope of the primordial power s pectrum n = 
0.93 - 0.031 ln(fc/fc ) with k = 0.05 Mpc" 1 dSpergel et alJ 
2003). Using these parameters we apply the rescaled ST ap- 
proximation to estimate the halo mass function in voids. 
Figure shows the estimates. As expected, the mass func- 
tion is then less steep than the mass functions with the same 
normalization and a constant slope n = 1. On the other 
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Figure 7. Comparison of the mass functions in five simulated 
voids (thin and thick solid lines correspond to voids with differ- 
ent radii as in Fig. |FJ to the predictions (dashed lines) based on 
constrained ST mass function (top panel) and rescaling the un- 
constrained ST mass function (bottom panel). The dashed curves 
from the highest to the lowest are for A(R vo ; ( j) = Cm, void = 
0.05,0.04,0.03. 



hand, decreased normalization makes it steeper so the net 
effect is that it is just shifted down with respect to the mass 
function found for as = 0.9 and n — 1. 

Note, that in Figure|5|the mass function is shown down 
to 10 h~ Mq. This is far beyond the range which can be 
tested in a numerical experiment, but due to the good agree- 
ment between the numerical simulation and the analytical 
approach we expect that the rescaled ST mass function can 
be extended to smaller masses. Using this approximation 
we predict that in a typical void of diameter ~ 20/i -1 Mpc 
there should be about 100 000 objects of mass greater than 
1O 7 /i _1 M0. It is a challenge to observers to confirm this 
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Figure 8. The solid curves are mass functions in simulations of 
voids. Void radii are reduced by 20% as compared with the void 
radii in Fig. [5] The two thick solid lines correspond here to the ra- 
dius -R V oid = 8/i _1 Mpc, while the three thin solid lines correspond 
to the radius -R vo jd = 6.4/i — ^^Mpc. Dashed curves show predic- 
tions from the rescaled ST mass functions for f^M.void = 0.03 and 
0.04. The analytical predictions and numerical results are now in 
a very good agreement. 
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Figure 9. Rescaled ST mass function for SIm vo j c i = 0.04 
and erg = 0.9 (dashed line) in comparison to the mass func- 
tion obtained for erg = 0.84 and scale-dependent spectral index 
n = 0.93 - 0.031 ln (fc/fcp) (solid line) as proposed by WMAP 
iSpergel et alj|2003f> . 
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Figure 10. The number of haloes with the virial mass larger than 
lO 12 fc, _1 M (bottom curve), 10 10 h,~ 1 M Q (middle curve), and 
10 8 /i — 1 Mq (top curve) inside a typical void of radius 10h~ 1 Mpc 
as a function of the average density of the void fi V oid . These pre- 
dictions come from the rescaled unconstrained mass function. A 
"void" with n vo i(j = 0.3 has the average matter density in the 
Universe. 

prediction. Recently, the existence of a lar ge population o f 
hydrogen clouds in voids has been claimed iMannindl2 002 ) . 
These clouds could be associated with a fraction of the nu- 
merous haloes which we found in the simulated voids. 

The expected number of dark matter haloes of different 
mass depends on the void (under) density. In Figure ITTil wc 
show predictions for the number of haloes inside a typical 
void of radius 10/i _1 Mpc for haloes with mass larger than 
10 12 h~ 1 M Q , lO lo /i _1 M and 1O s /i _1 M as functions of the 
average density of the void. There is a dramatic decline in 
the number of haloes when the average density of the void 
falls below a certain percentage of the mean density. The 
larger the mass of the haloes, the larger is the void density 
at which the sharp decline happens. This decline is related 
to the steepening of the void mass function at the high mass 
end as seen in Figure |S| Note, that for typical voids with 
the average density ten times below the mean density, the 
number of 1O 8 /i _1 M haloes is not yet suppressed: the sharp 
decline is at lower average densities. 



6 DISCUSSION AND CONCLUSIONS 

We used iV-body simulations to study the formation of voids 
in the large-scale structure of the Universe. We first identi- 
fied dark matter haloes and then searched for voids in the 
distribution of the haloes. We found that the size of the 
voids depends weakly on the lower mass (or circular veloc- 
ity) limit of the haloes chosen to define the void. Voids in 
a set of more massive haloes tend to be slightly larger than 
voids in a set of smaller haloes. 

The interior of five largest voids was resimulated with 



very high mass resolution of 4x 1O 7 /i -1 M . These voids have 
diameters of about 20/i -1 Mpc and their density is about 
a factor of 10 smaller than the mean density in the Uni- 
verse. Inside a void of this size we found typically about 
50 haloes with circular velocities larger than 50 km/s and 
more than 800 with circular velocities larger than 20 km/s. 
A scale-dependent slope of the primordial p ower spectrum 
as recently suggested bv lSpergel et alJ (|2003j) would slightly 
red uce the number of low m ass haloes in voids. 

iMathis fc White! <|2002|) find "several void regions with 
diameter 10/i _1 Mpc in the simulations where gravity seems 
to have swept away even the smallest haloes" they "were 
able to track". According to their Table 1 they tracked 
haloes down to 10 particles corresponding to halo masses 
3.6 x 10 10 /i- 1 M Q . In our simulations of voids of diame- 
ter 20h~ Mpc we find typically up to 10 haloes of this 
mass (represented by 1000 particles) . More than half of the 
haloes are close to the outer void boundary with a dis- 
tance to the void's centre larger than 80% of the void's 
radius. Since our void volume is 8 times larger and more 
massive haloes tend to be situated in the outer part of the 
void it would be not extremely unprobable to find inside 
our voids regions of 10h~ Mpc free of any halo more mas- 
sive than 3.6 x 10 10 h~ 1 Mp,. In this respect we agree with 
IMathis fc White! (I2002T) . Following the evolution of voids nu- 
merically it seems that these haloes are not swept away by 
gravity but never form in regions with the lowest density. 

The fact that regions devoid of haloes with masses 
larger than 5 x lO lo /i _1 M are smaller than those devoid 
of haloes with masses bigger than 2 x 10 11 h~ 1 M empha- 
sizes our argument that the size of the voids depends on 
the objects used to define the void. For example, the voids 
in the distribution of clusters are not empty: they contain 
many L* galaxies. Likewise, the voids in a sample of L„ 
galaxies should contain dwarfs. This argument suggests that 
voids in the dark matter distribution should be self-similar 
in the same sense as clusters of galaxies have hundreds of 
galaxies and galaxies have hundreds of satellites. It must be 
so as long as only the gravity is the dominant factor and 
the spectrum of fluctuations is approximately scale invari- 
ant. The question whether luminous galaxies form or do not 
form in all these dark matter haloes is then the most impor- 
tant question in the theory of galaxy formation in different 
environments. 

Assuming with Mat his fc White! (Hp02) a luminosity 
Mb = —16.5 for a galaxy hosted by a halo of 3.6 x 
lO lo /i _1 M we predict about five of these galaxies to be 
found in the inner part of a typical void of diameter 
20/i _1 Mpc. In principle they can b e detected. In t he gi- 
ant (radius 31.5 Mpc) Bootes void ISzomoru et alJ (Il996f) 
study void galaxies in about 1% of the void volume. All but 
one of them were substantially brighter than M = —16.5. 
The one reported with M — —16.2 is the companion of a 
brighter one. To compare the model predictions with obser- 
vations one would have to study a void in the distribution 
of galaxies with limiting magnitude Mb between —17.5 and 
— 18.5 which roughly corresponds to our threshold mass of 
2 x 1O 11 /i~ 1 M . 

Note, that so far we assumed that each dark matter 
halo hosts a galaxy. This may not be true. Physical processes 
of galaxy formation are not well known and there could be 
processes that strongly suppress the formation of stars inside 
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small haloes, which collapse relatively late in voids. One such 
proc ess which has been widely discussed is t he ionizing flux 
(e.g. iBullock et afll200ri : iBenson et afll2002h . 

We determined the mass function of dark matter haloes 
in the simulated voids and compared it with the ana- 
lytical predictions based on the Sheth-Tormen formalism. 
The formalism was applied in two versions: (1) We used 
the ansatz for the con strained mass function proposed by 
ISheth fc Tormenl l|2002h and (2) we propos ed our own exten- 
sion o f the unconstrained mass function of ISheth fc Tormenl 
( 1999) by rescaling the power spectrum in the void. We 
found that our approach is not only much simpler in applica- 
tion, but also reproduces more accurately the mass functions 
obtained in the simulations of voids. 
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